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ABSTRACT 



We explore the enigmatic population of long-period , apparently non-rec ycled 
pulsars in globular clusters, building on recent work by iBoyles et al.l (120111 ) This 
population is difficult to explain if it formed through typical core collapse su- 
pernovae, lea ding many authors to invoke electron capture supernovae. Where 
Boyles et al.l dealt only with non-recycled pulsars in clusters, we focus on the pul- 
sars that originated in clusters but then escaped into the field of the Galaxy due 
to the kicks they receive at birth. The magnitude of the kick induced by electron 
capture supernovae is not well known, so we explore various models for the kick 
velocity distribution and size of the population. The most realistic models are 
those where the kick velocity is < 10 kms -1 and where the number of pulsars 
scales with the luminosity of the cluster (as a proxy for cluster mass). This is 
in good agreement with other estimates of the electron capture supernovae kick 
velocity. We simulate a number of large-area pulsar surveys to determine if a 
population of pulsars originating in clusters could be identified as being separate 
from normal disk pulsars. We find that the spatial and kinematical properties of 
the population could be used, but only if large numbers of pulsars are detected. 
In fact, even the most optimistic surveys carried out with the future Square Kilo- 
meter Array are likely to detect < 10% of the total population, so the prospects 
for identifying these as a separate group of pulsars are presently poor. 
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Introduction 



There are currently 144 pulsars^ known in 28 Galactic globular clusters (GCs). The 
vast majority of these are millisecond pulsars (MSPs), characterized by short spin periods 
and small period-derivatives. Such a population arises quite naturally in GCs, which are 
old stellar system s that are thought to contain a reservoir of primordial neutron stars (NSs) 
(IHut et al.lll992l ). In the dense environments in the cores of most GCs, these NSs undergo 
frequent exchange interact ions, and may the n be "recycled" into MSPs by accreting matter 
from a binary companion (lAlpar et al.lll982l ). In addition to these MSPs, however, there is 
a small but enigmatic population of long-period pulsars that seem similar to the "normal" 



pulsars usually found in the Galac tic disk (see Figured]) ( iBiggs et al.lll994j ; iLyne et al.lll996 



Chandler 


2003: 


Lvnch et al. 


2012) 



throughout this paper to distinguish them from the normal Galactic disk pulsars. The 
standard scenario for forming normal disk pulsars involves the core collapse of a massive 
star, giving rise to a young pulsar with a h igh magnetic field, which quickly spins down to 

The inferred lifetimes of normal pulsars are 



P ~ 0.3 s ( Faucher-Giguere fc KaspJ 



typically 10 to 100 Myr ( iRidley fc Lorimer 



2006 



2010J), and as such, they are usually associated 



with the recent death of massive stars, which themselves have lifetimes < 100 Myr. This is 
where the mystery of NRPs in GCs arises — GCs are composed of old, < 1 M stars, and all 
stars massive enough to form NSs (along with any pulsars that were formed) should have 
died some ~ 10 Gyr ago. The fact that several NRPs are observed in GCs requires an 
alternative to the standard core collapse model. 



Several authors ( ILyne et al.l Il996l ; llvanova et al.l 120081 ) have invoked the collapse of a 
massive O-Ne- Mg white dwarf (WD ) via electron capture (so called electron capture super- 
novae, or ECS (lNomotolll984l . 119871 )). Recent theoretical work supports the no tion that ECS 



are e ssential for understanding the full population of neutron stars in GCs ( llvanova et al. 



20081 ). Unfortunately, ECS have never been directly observed and their energetics remain 
uncertain, although there is good reason to believe that they are about an order of magni 



tude less energetic than core collapse supernovae (CCS) (IKitaura et al.ll2006l ; iDessart et al. 
20061 ). This probably leads to sm all natal kicks when com bined with a faster and more 



symmetric explosion than in CCS (IPodsiadlowski et al.l 120041 ) . NRPs may be an important 
observational constraint on the properties of ECS. 



Motivated by the lar ge number of recent and sensitive pulsar searches of GCs, iBoyles et al. 



( 1201 ll . hereafter iPaper J ) used statistical techniques to explore the underlying population of 
NRPs in GCs. The results are based solely on observations and some simplifying assump- 



4 For an up-to-date list see http://www.naic.edu/~pfreire/GCpsr.html 
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tions about the luminosity function and lifetime of NRPs, and as such provide a constraint 
on the birthrate of NRPs in GCs, regardless of how they are formed. One of the key results 
from this study was the compilation of probability distributions for the birthrates of NRPs 
in the majority of GCs. 



Paper II gave detailed consideration only to those pulsars that gained a sufficiently small 



natal kick as to be retained by their host clusters. Given the high upper limits on the 
birthrate of GC NRPs and the relatively shallow potentials of most GCs, there is an intriguing 
possibility that a large population of NRPs may have escaped from their progenitor clusters 
at birth and entered the field of the Galaxy, where they coul d contri bute to the observed 
population of normal disk pulsars. Building on the results of iPaper J . we have carried out 
Monte Carlo simulations to explore the properties of this purported population, and the 
feasibility of detecting it as a separa te group. In §2] we provide a brief overview of the 
techniques and results from IPaper t . In £J3] we describe our simulations and present the 
results in £|H We discuss the implications of these results in £J5] and summarize in £j6j 



2. Overview of IPaper I 



The goal of IPaper t was to characterize the underlying population of NRPs in a cluster 
using observational constr aints. Ob viously, this approach is only valid for clusters that have 
been sea rched for pul sars. IPaper II compiled results from se arches of 97 clusters (out of 156 
listed in lHarrisI (120101 )) carried out by numerous groups (see lPaper II for a complete list and 
references). For each search, the detection probability was defined as 



9(L min ,f(L)) 



X 



where L min is the limiting luminosity of the search, and f(L) is the luminosity distribution of 
NRPs. In words, 9 is simply the fraction of pulsars that lie above L m\ n . Five log- normal distri- 
bution s with differe nt parameters were used for f(L); one given by lFaucher-Giguere fc Kaspi 
( 120061 ) and four by iRidley fc Lorimerl ( 120101 ). This gave rise to a range of values for 9. 



Bayes' theorem was then used to obtain a probability density function (PDF) for the 
total number of potentially observable NRPs (N), given the number of actual detections 
(n), and 9. Potentially observable means all NRPs that currently reside in clusters and 
whose radio beams cross our line of sight (i.e., all those that could be observed with infinite 
sensitivity). We shall discuss corrections to this number in §21 This PDF is 



V{N, 9\n) oc C(n, N\6)V(N, 9), 



(2) 
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where C(n, N\9) is the likelihood of detecting n pulsars from a population of N given 9, and 
V(N, 9) is the joint prior PD F for N and 9. Since pulsar searches have only two possible 
outcomes (success or failure), iPaper t chose the binomial distribution for the likelihood: 

£( "^ |S) = S(F^)! S " (1 - S) "" (3) 



To characterize the joint prior, V(N,9), IPaper II first assumed that N and 9 were inde- 
pendent. The prior for N was taken to be a uniform distribution on [n, oo), so that the 
only restriction is N > n. The prior for 9 was also taken to be uniformly distributed on 
[9 m m , 9 max \, where ^min and # max correspond to the lowe st and highest values of 9 for different 
choices of f{L). Having obtained V(N, 9\n), IPaper t then marginalized over 9 to obtain a 
final PDF for N (denoted as V(N\n)) for each of the 97 clusters in the sample. 



3. Simulating NRPs that Escape from Clusters 

As mentioned in £j2j it is necessary to correct V(N\n) to account for the fraction of 
NRPs whose radio beams do not cross our line of sight, and for those who are no longer 
bound to their progenitor GCs. This second correction is necessary because GCs have fairly 
shallow potentials, with escape velocities t> esc ~ 30 kins" 1 . Some NRPs will inevitably 
escape the cluster due to the natal kicks they receive during formation and enter the field of 
the Galaxy. This is precisely the population that we are interested in here. We have used 
Monte Carlo simulations to model the evolution of NRPs that escaped from GCs and to 
determine if this population could be distinguished observationally from other pulsars in the 
Galaxy. We proceed in four steps: i) we determine the number of NRPs that will enter the 
Galaxy; ii) we model the spatial and kinematical evolution of this population; Hi) we model 
the rotational and electromagnetic evolution; and iv) we "detect" this evolved population 
by simulating a number of large-area pulsar surveys. As described in the following sections, 
we actually explore several different types of models for the population, and for each model 
we also explore several natal kick distributions. To ensure robust statistical results, each 
combination of model class and kick distribution is simulated ten times, and the final results 
are taken to be the mean of these runs. Figure [2] is a schematic overview of our simulations. 
We discuss each step below. 



3.1. Determining the Number of Escaped NRPs 



We choose to characterize V(N\n) by the upper bound of the 95% confidence interval, 
which we denote as N 95 (note that the birth rates discussed in IPaper j used the median of 
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this distribution). This is corrected for beaming and escaping pulsars, and divided by the 
mean lifetime of NRPs to obtain a birth rate for each cluster that we consider, 



K 



or, 



VbeamVcTv ^NRP 



(4) 



where 1Z„ V is the birth rate, ^beam = 0.1 is the beaming fraction ( jTauris fc Manchesterlll9981 ). 
r) av is the correction for escaping pulsars, and r NRP is the average pulsar lifetime. We derive 
an average lifetime of 43 Myr by taking the total number of radi o-loud pulsars (1.2 x 10 6 ) 
and d ividing it by the Galactic pulsar birthrate (2.8 century -1 ) ( iFaucher-Giguere fc Kaspi 
20061 ). The subscript a v refers to the mean birth velocity dispersion. The actual value of 
o v for NRPs is unknown, since they are probably not formed via typical CCSNe , so we 
vary this parameter in our simulations. As we shall see, there is good reason to believe 
a v ~ 10 kins" 1 . We assume that birth velocities can be described by a Maxwell-Boltzmann 
distribution so that the fraction of NRPs with v < v csc is 



fo° sc /mb(u) dv 
Jo^/mbM dv 

. , Vpxc i / 2 Dp 

erf 



7T (T T 



exp 



2^ 2 



(5) 



where /mb(^) is the Maxwell-Boltzmann distribution and erf signifies the error function 



erf (x) 



"* dt 



(6) 



Escape velocities for each cluster were taken from lGnedin et al.l (120021 ) . For clusters without 
a reported va l ue, w e assumed t> e sc — 30 kms -1 , which is roughly the median value from 
Gnedin et all J2OO2L 



Having obtained lZ av , all that remains is to choose a time scale, r max , over which we 
will consider the evolution of this population. The number of pulsars to simulate per cluster 
then becomes iV sim = 7Z av x r max . In our simulations, we choose r max = 200 Myr. As we 
explain in §3.3[ many pulsars will cease radio emission on timescales shorter than this. We 
choose a large value for r max to ensure that we treat long-lived NRPs properly. 

We refer to the above a pproach a s our "standard model" (hereafter Sa). A^ sim is obtained 
directly from the results of iPaper J . which is to say, without using any information about 
the clusters except for their escape velocities. However, the value of N 95 , and hence A^ sim , 
depends heavily on L m i n , in the sense that N95 will be larger (and not very constraining) for 
shallowly searched GCs. Therefore, we included three refinements to Sa in our study. 
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The first we call the "modified standard model" (Sb). For this, we simply exclude 16 
clusters with L 14 GHz,min > 7.5 mJykpc 2 , which reduces the number of shallowly searched 
clusters with very high values of N 95 . While this sensitivity cut-off is somewhat arbitrary, it 
strikes a good balance between keeping enough clusters for the simulations to have meaningful 
results, while preventing unwieldy computations due to extremely large numbers of simulated 
escaping NRPs. 

In the next two models, we choose a reference cluster, and scale all other values of N sim 
by some chosen parameter. We use M22 as a reference cluster because it has the lowest 
value of L min and consequently the most tightly constrained N S i m (relevant parameters of 
M22 can be found in Tabled]). In the first approach, we scale by the V-band luminosity of 
each cluster, L v , so that 

^shn = -fysim,M22 J ^ ■ (7) 
-^V,M22 

This is the "luminosity scaling model" (L). Our reasoning is that at the very least, the 
number of NRPs produced in a cluster should scale with the number of progenitors in the 
cluster. However, to determine the number of progenitors would require first identifying 
those progenitors, and then making some assumptions about how their number scales with 
other properties of the cluster. We sought a more general approach, hence our decision to 
use the luminosity of the cluster — brighter clusters should be more massive (to within a 
factor of the cluster's mass-to-light ratio), and more massive clusters should contain more 
NRP progenitors. Furthermore, the V-band luminosity is fairly easily determined as long as 
the distance to the cluster and reddening effects are well constrained!^, and is also readily 
available for nearly all Milky Way GCs. 

Finally, we also scale by the so-called core interaction rate, T c oc pl' 5 r^, where p c is the 
central density of the cluster and r c is the core radius; thus we have 

r c 

N sim = -Nsim,M22=; . (8) 

1 c,M22 

We call these the "interaction rate scaling models" (G). T c has been shown to correlate 



well with the num ber of low- mass X-ray binaries and MSPs in a cluster (jPooley et al.l 12003 



Abdo et al.ll2010l ). further supporting the notion that both populations are related to binary 
exchange interactions. A leading explanation for NRPs in GCs are ECS, particularly via 
accretion or merger induced collapse. These scenarios also require mass-transfer or merger 



5 The values we use are from 



Harris 



(|2010h and have been corrected for reddening. 
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in a binary system. Hence, it is interesting to explore models in which NRPs also have some 
dependence on r c . 

For each of these models, relevant parameters (iV95, Ly, and T c ) were only available 
for a subset of all 156 Milky Way GCs. It will be necessary to adjust our final results to 
account for unmodeled clusters, but we save this step for §U Table |2] summarizes each class 
of models. 



3.2. Spatial and Kinematical Evolution 

In this step, we model the evolution of the escaping NRPs as they travel through the 
Galaxy. We begin by defining a Galactocentric coordinate system with 

x = D cosbsmi, (9) 
y = D Q — D cos b cos £, and (10) 
z = Dsinb (11) 

where D is the distance, D Q = 8.5 kpc is the Sun's distance from the Galactic center, b 
is Galactic latitude, and £ is Galactic longitude. Note that this differs from the typical 
Cartesian Galactic coordinates, where it is usually the x-axis that connects the Sun and 
Galactic center. This was done for compatibility with previously developed software tools. 

We assign each NRP initial coordinates that are equal to the coordinates of the host 
GC. Initial 3-D velocity components (v x ^, v y ^ v z ^) are chosen at random from a normal 
distribution with a zero mean and standard deviation of a v (this is equivalent to a Maxwell- 
Boltzmann distribution for the full space velocity). If the total space velocity with respect 
to the GC, v\ = (v^i + Vy { + v ^ ; ) 1//2 > v esc , then we consider the pulsar to have escaped from 
its host cluster and entered the field of the Galaxy. Although we keep track of the number 
of pulsars retained by their host clusters, we do not consider them further. All subsequent 
analysis on the number of emiting and detectable NRPs deals only with those that escape. 
We treat each GC as being fixed at its current position in the Galaxy. In reality, the clusters 
are on their own orbits, and the pulsars will inherit the systemic velocity of their progenitor 
clusters. Each pulsars enters the field of the Galaxy with Vfi eld = (vf — v^) 1 ^ 2 . 



Having chosen initial spatial and velocity c omponents, we integrate th e m otion of each 

pulsar through the Galactic potential. Following lCarlberg &: Innanenl (119871 ) and lKuijken &: Gilmore 
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( 119891 ) (see also lWex et al.ll2000l ). we use a 3-component Galactic potential of the form 



l (R,Z) 



-GM 1 



+ b l + R 2 



-1/2 



;i2) 



where the superscript i = D, B, N indicates the contribution from the disk-halo, bulge, or 
nucleus, respectively. The full potential is the sum of these three components. The values 
of M\ a 1 , b\ (3j, and hj can be found in Table [3) The integration time, t em \, is chosen at 
random from a uniform distribution on the interval [0, 200 Myr]. In other words, we assume 
that an NRP may be born at any point in the past 200 Myr, and evolve it forward to the 
present. This also assumes that the birthrate of NRPs is constant over this interval. The 
final spatial and kinematical properties are then recorded for later analysis. 



3.3. Rotational Evolution and Energetics 



Birth spin periods and surface magnetic fields (BS ) are chosen at random for each 
escaping NRP. Following IFaucher-Giguere fe Kaspil (120061 ) . we use a normal distribution for 
P with a mean of 0.3 s and standard deviation of 0.15 s. Negative periods are rejected and 
redrawn. We also use a log-normal distribution for B s , with a mean in the base- 10 logarithm 
of 12.65 and standard deviation 0.55. 

We evolve the pulsar's rotation under the assumption of pure magnetic dipole braking 
and a constant magnetic field, so that the observed period and period derivative are 



P 2 + lQTtB 2 RH evol 



3c 3 / 



1/2 



and 



P 2 - P 2 



(13) 
(14) 



2-P^evol 

where R = 10 km and / = 10 45 gmcm 2 are the assumed radius and moment of inertia for 
the pulsar, respectively, and t evo i is defined as above. Hence, we arrive at the final P and P 
at the end of our simulation. 

We calculate the observed luminosity using a power-law model that depends on P and 
P. We prefer this model because it relates the luminosity to the ro tational energy loss of 
the pulsar. Once again, we turn to IFaucher-Giguere &: Kaspil (120061 ) for the exact form of 
this power-law: 

1/2 



L = 10 icor L P-' S/2 



P 



lO" 15 



ss 



(15) 
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Lq = 0.18 mJykpc 2 is a "standard" luminosity and L cor is a correction factor that accounts 
for uncertainty in the model. It is drawn at random from a normal distribution with zero 
mean and standard deviation <Jl coi — 0.8. 

Pulsars have finite lifetimes, and at some point will cease radio emission. This seems 
to correspond to a "death -line" in the P-P diagram, which is well described theoretically 



(IBhattacharya et al.lll992l ) and empirically as 

|| = 1.7xlO n Gs- 2 . (16) 
pi 

Any pulsars with B S P~ 2 less than this value will no longer be visible. We track the evolution 
of these pulsars, but do not include them for consideration in the next step. 



3.4. Simulated Surveys 



The final step is to "search" for potentially visible NRPs by simulating various large- 
area surveys. A pulsar can be detected only if its sig nal-to- noise ratio (S/N) li es above the 
minimum S/N threshold of a given survey. Following lLorimer fc Kramer! (120041 ). 



S/N = S, 



/3(T sys + T E 



sky J 



P 

w- 1 



(17) 



where S u is the flux density of the pulsar, G is the telescope gain, iVp i is the number 
of summed polarizations, Av is the bandwidth, t int is the integration time, (3 is a factor 
that accounts for quantization losses, T sys and T s k y are the system and sky temperatures, 
respectively, and W is the observed pulse width. Most of these factors are intrinsic to the 
specific survey in question, but T sky and W depend on the position and properties of the 
p ulsar. Sky temperatu res at the position of each pulsar are taken from the 408 MHz survey 
of lHaslam et al.l ( 119821 ) and scaled to the appropriate frequency assuming a power-law with 
a spectral index of —2.6. The pulse width is described as the quadrature sum of 



W 2 = W 2 + + t 2 

' ' " ' int 1 samp 



where W- m t is the intrinsic pulse width, t S am P is the instrumental sampling time, £dm is the 
dispersive smearing within a given frequency channel, and r s is the scattering time. We 
model the intrinsic width as 



log W mt = log 



0.06 



P 



ms 



0.9 



(19) 
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whe re q is a random vari able drawn from a normal distribution with a standard deviation of 
0.3 (ILorimer et al.l 120061 ) . The sampling time for each survey is listed in Table HI The dis- 
persion measure (DM) is calculated usin g the known distance t o the pulsar and the NE2001 
model of Galactic free electron density (jCordes &: Laziol 120021 ) . The dispersive smearing is 
then simply 



"DM 



8.3 x 10 3 s 



DM 



pc cm" 



Au, 



chan 



MHz 



v 



MHz 



(20) 



where Az/ C han is the channel width and v is the center frequency ( ILorimer fc Kramer! 120041 ) . 
Various empirical relations betw een r s and DM appear in the literature. We use the rela- 
tionship given by ICordesl ( 120021 ) : 



log — 



Ts 
flS 



-3.59 + 0.129 log 



DM 



pc cm 



+ 1.02 



log 



DM 



pc cm 



-3 



4.4 log 



v 



GHz 



(21) 



One factor that we do not model here is the affect of radio frequency interference (RFI), 
which can be particularly problematic for blind searches of long-period pulsars. Obviously, 
some surveys will be more affected by RFI because of proximity to man-made sources and/or 
poor instrumental resistance to strong RFI. Nonetheless, careful data analysis can help to 
mitigate these effects. 

After determining the relevant parameters for each NRP, we use the survey tool from the 
psrpop software packagtfl (ILorimer et allfoiofih to simulate the surveys. Pulsar luminosities 
are scaled to the appropriate frequency assuming a power-law spectral index of —1.6. 



We simulate five surv eys: the Parkes M ultibeam Survey (PMSURV) (ILorimer et al. 



20061 ): the P-ALFA survey (ICordes et all 120061 ). the GBT North Celestial Cap Survey (GB- 
NCC); a visible-sky GBT survey similar to the GBNCC survey (GB TALL); and a hyp othet- 
ical all-sky survey using a future Square Kilometer Array (SKA) (jSmits et al. 20091 ). The 
characteristics of each survey can be found in Table HI Each survey is simulated 100 times, 
for each run of the simulation, allowing us to characterize the median number of detected 
pulsars. Because the spatial, kinematical, and rotational evolution of the NRPs is simulated 
ten times for each combination of model class and a v (see above), we have ten values for the 
median number of detected pulsars from each combination. The final reported number of 
detected NRPs (iV det ) is the mean of these ten values plus the standard error, and represents 
an upper limit. 



e http : //psrpop . sourcef orge . net/ 
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4. 



Results 



The results of our simulations can be found in Table |5j The values reported here have 
been multiplied by a scale factor (156/A^gc) that accounts for clusters not included in that 
model class. This is equivalent to assuming that un-modeled clusters follow the average 
results for that model class. The total number of all NRPs that escape from their progenitor 
clusters (N esc ) and that are retained by their host clusters (iV re t) are listed in the first two 
columns. These include pulsars whose radio beams do not cross our line of sight as well as 
those that have crossed the death line; these will be dealt with in §4.21 Hence, these numbers 
represent 95% confidence interval upper limits on the total number of NRPs formed in the 



It is immediately obvious that, for a given model class, iV esc increases with increasing 
a v . This is not at all surprising, since a larger kick velocity dispersion will lead to more 
high velocity pulsars that can escape the cluster potential. In contrast, N Tet remains nearly 
constant with o v . This is because we have essentially normalized all values by the "observ- 
able" (i.e., based on observations) value of Ng§, which characterizes the number of retained 
NRPs. We point out that 7Z av also increases with a v , so that there are more pulsars in 
total simulated for higher velocity dispersions. Model class Sa gives rise to the largest total 
number of NRPs by far. This is entirely due to the very high values of N 95 for most GCs. 
Even when we exclude the most unconstraining clusters from consideration in model class 
Sb, the total number of NRPs is still very large. Model class G produces fewer pulsars than 
Sb, while model class L produces the fewest of all. As we shall see in §5.1[ this has important 
consequences for the viability of ECS as an explanation for NRPs. For now, we turn our 
attention to exploring the characteristics of the escaped population in more detail. 



Figure [3] shows a sample sky distribution of all escaping NRPs for model L70 in Galactic 
coordinates (other models have very similar forms with different numbers of pulsars). The 
pulsars tend to group around their progenitor clusters. In models with higher a v the pulsars 
travel further from their host clusters, as expected, but there are also more pulsars in total. 

We also model calculate the proper motion of each pulsar. The proper motions are a 
combination of the kick velocity of the pulsar (evolved as the pulsar travles through the 



N esc and iV re t scale linearly with the choice of r max , so it is trivial to adjust these when considering a 
different timescale. 




4.1. Spatial and Kinematical Properties 
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Galaxy) and the systemic velocity of the GC itself (both projected onto the plane of the 
sky). Pulsars from a common progenitor cluster will thus have similar systemic velocity 
components. The velocities of many GCs arepdf observed to be ~ 100 kms -1 , larger than 
the kick velocities received by the NRPs at birth. 

4.2. Number of Detectable NRPs in Surveys 

The results of our simulated surveys can be found in Table [51 It is at this stage that 
we account for pulsars that have crossed the death line. The third column gives the number 
of all radio emitting pulsars (iV emit ). The number of detected pulsars for various surveys is 
given in the following column^. These numbers have already been multiplied by a constant 
Vbeam = 0.1 to account for pulsars whose radio beams do not cross our line of sight. It is 
immediately obvious that significant numbers of GC NRPs are detected only when N em n 
is large (corresponding to large a v ). The most successful survey is, not surprisingly, our 
hypothetical all-sky SKA survey. Such a survey would presumably be designed to detect 
nearly all potentially observable disk pulsars in the Galaxy, but only ~ 8% of emitting GC 
NRPs are detected. The next most promising survey is the hypothetical GBTALL, but this 
detects only ~ 0.03% of emitting NRPs. The success of the PMSURV and P-ALFA are less 
than but comparable to the GBTALL. The GBNCC detects very few, if any, NRPs. We can 
attribute this to the design of the survey, whose goal it is to find nearby and bright MSPs; 
the 1374 MHz limiting flux density is an order of magnitude higher than PMSURV, P-ALFA, 
or the SKA survey. However, the GBTALL has an identical set-up and does much better. 
We attribute this to the large survey area and better coverage at low declinations, where 
most GCs are found. The PMSURV and P-ALFA both utilize high observing frequencies 
and long integration times to search for highly dispersed, and thus generally more distant 
pulsars. At first glance it may seem strange that so few NRPs are detected, given that we 
already know of four that reside in GCs. These pulsars were detected in targeted surveys 
with very long integration times. Furthermore, for low a v , we predict N esc < N Tet ; in models 
where N esc w N Tet , the number of detected NRPs is closer to what is actually observed in 
clusters. 

We can conclude that, to have any hope of detecting large numbers of escaped NRPs, 
surveys must be very sensitive and cover as much sky as possible. Even then, success will 
depend strongly on the size of the underlying population, which depends on the kick velocity 



s N em n and Ndot are not very sensitive to the choice of r ma x unless it is less than the lifetime of a typical 
NRP. 
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of NRPs. 



5. Discussion 
5.1. Which Models are Realistic? 



One of the primary findings of iPaper II was that very high values of a v greatly over- 
produce NRPs in GCs according to any reasonable metric. For example, they find a birth 
rate of 422 and 0.25 psr century" 1 for a v = 130 and 10 kms -1 , respectively (keep in mind 
that these are for the median number of NRPs, whereas we use the higher 95% confidence 
upper limit). For comparison, the birth rate of normal disk pulsars is 2.8 psr century" 1 



(IFaucher-Giguere &: Kaspil 120061 ) . More relevant for the case of NRPs is how these results 
compare with the population of WDs that are the likely progenitors to NRPs via ECS. 
Paper I estimated the total number of WDs that could potentially form NSs via ECS0 to be 



~ 2.1 x 10 5 across the entire GC population. If we assume that all of these WDs eventually 
become NSs and that the birthrate of NRPs is constant, then the timescale for exhausting 
this population is simply 2.4 x lQ b /lZ av . These timescales, r exhaust , are given for each of our 
models in Table Even for the smallest lZ av in our standard models, T ex h aU st ^100 Myr. 
This is only ~ 1% the age of a typical GC (~ 10 Gyr). Although these birth rates are 
only upper limits, such a large discrepancy seems difficult to explain. Recall, however, that 
our standard models are heavily affected by the limiting luminosity of cluster searches, so 
we believe that these results simply indicate that the standard models are not very well 
constrained. Model class G, where we scale the number of NRPs formed in M22 by the core 
interaction rate of a cluster, produces fewer NRPs than our standard models. Nonetheless, 
inspection of Table [6] indicates that the implied birthrates are still probably too high. 

However, the situation improves for model class L, where we scale by the luminosity 
of a cluster. High velocity dispersion cases still seem to overproduce the number of NRPs 
considerably, but for a v < 30 kms" 1 the implied birth rates are low enough that the resevoir 
of WDs may last for > 15% of the age of the typical cluster. When one considers that the 
birth rates reported here are only 95% confidence upper limits, it seems that ECS are a viable 
explanation for NRPs after all. Given that the true efficiency in going from WDs to NRPs is 
probably < 100%, we favor a typi cal ECS kick veloc ity of < 10 kms" 1 . Other authors have 



inferred small o v for ECS as well. iPfahl et al.l (120021 ) studied the eccentricities of high-mass 



X-ray binaries that are believed to form via ECS and concluded that o v < 50 kms x . More 



3 Assumed to be all WDs with 1.0 < M < 1.4 M 
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recently, iMartin et all (120091 ) used a similar line of reasoning to argue for a v ~ 15 kins" 1 
(though the authors note that this may be in conflict with the observed misalignment between 
the rotational and orbital spin axes in these systems ). PS R J0737— 30 3 9B is also believed 
to have formed via ECS JPodsiadlowski et all l2005h and IWong et all feoioh find a 95% 
confidence interval of 5-120 kms -1 for the kick of this neutron star. Our results are in 
general agreement with the notion that ECS must be less energetic and lead to smaller 
neutron star kicks than CCS, but may be more constraining than other estimates. 

With a firm estimate of the efficiency for converting WDs to NRPs and the duration 
of this process, a more exact estimate could be made. There are two important caveats to 
keep in mind, however. The first is that model classes L and G are based on the upper limit 
of the number of NRPs in only one cluster, M22. We chose this cluster because it has been 
searched to a lower limiting luminosity than any other. If M22 is an atypical cluster in terms 
of its NRP content, then this would bias our results. However, other clusters (e.g. M28, 
47 Tucanae, Terzan 5) have been deeply searched and also have small implied NRP birth 
rates. Model classes L and G also a ssume th at only one characteristic of a GC influences the 
number of NRPs that are formed. iPaper t found evidence suggesting that metallicity may 
play a significant role in the formation of NRPs, and we cannot rule out other factors. 



The best way to differentiate between these results and those presented in lPaper II would 
be to search a large nu mber of G Cs more deeply. This would yield much tighter constraints 
for the method used in IPaper J . and would reduce the number of simplifying assumptions 
we need to make here. 



5.2. Can the Population of Escaped NRPs Be Identified? 

The next question to ask is, what are the prospects for identifying the population of 
escaped NRPs as separate from normal disk pulsars? In the following discussion, we will limit 
ourselves to the results for model class L. With large numbers of pulsars it may be able to 
identify NRPs based on their spatial distribution around GCs or their kinematic similarities, 
since NRPs born from the same GC should have the same systemic contributions to their 
velocities. Unfortunately, we favor models with low a v , and as §4.21 makes clear, in this case 
even the best surveys can detect < 100 NRPs. Proper motions would probably only be 
available on a sub-set of these pulsars. However, when the SKA comes online, it will be able 
to measure proper motions interferometricaly on nearly all pulsars it detects, so this may be 
a future avenue for obtaining proper motions. Nonetheless, it is worth keeping our results 
in mind as high sensitivity surveys are carried out. 
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Large area surveys are not the only means by which GC NRPs could be detected. Since 
many escaped NRPs are found to remain fairly close to their progenitor clusters, especially 
for our favored low a v models, targeted surveys of the regions around GCs could be more 
fruitful, as it would allow for deeper integrations and better sensitivity We have calculated 
the minimum distance between an escaped NRP and a GC for all emitting pulsars in our 
simulation. There is some dependence on <r v , but in general ~ 3-7%, ~ 13-16%, and ~ 20- 
25% of NRPs lie within 1°, 3°, and 5° of a GC, respectively. However, a large a v would 
still give rise to a larger pool of potentially observable pulsars. A survey covering a ~ 3° 
radius around all GCs would then, in principle, capture roughly 15% of the population of 
NRPs. However, with 156 Milky Way GCs, this would amounts to a total survey area of 
over 4400 square degrees, with each portion requiring substantial integration times. Such a 
survey would be very difficult. The situation is improved if such a survey were limited to 
only 1° around each cluster (reducing the survey area to about 500 square degrees) or by 
targeting only the most promising GCs, such as those already known to contain an NRP, 
but this would also decrease the likely number of NRPs that could be detected. 



Conclusions 



Paper II explored the enigmatic population of NRPs found in GCs by setting limits on 
the size of the population using the results of various GC surveys. We have built upon these 
results to explore the prope rties of N RPs that escape from clusters and enter the field of 



the Galaxy. We agree with iPaper II that current GC surveys are not sensitive enough to 
constrain the population on their own. However, if we assume that the number of NRPs in 
a GC scale with some properties of the cluster, specifically luminosity (as a proxy for cluster 
mass), then the size of the population is more realistic. It seems that in these cases, ECS 
may be sufficient to account for GC NRPs. We favor models that rely on luminosity scaling 
and invoke a low natal kick velocity (~ 10 kms -1 ), but there is still too much uncertainty in 
how NRPs are formed to place more definite limits. Unfortunately, large numbers of escaped 
NRPs would be difficult to detect with any current large-area pulsar surveys, so the chances 
of identifying them as a separate population are slim. The best prospects lie with future 
surveys with the SKA. If sufficiently large numbers of GC NRPs can be detected in the field 
of the Galaxy, it may be possible to distinguish them from normal disk pulsars through their 
spatial distribution and proper motions. The best approach for learning more about NRPs, 
though, remains identifying those that are bound to their progenitor GCs. This will require 
more sensitive searches of many clusters. 



We are extremely grateful to Norbert Wex for providing us with an independent method 
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Period (s) 



Fig. 1. — The P-P diagram showing the four NRPs in GCs as b lack triangles. The gra y 
points are all Galactic pulsars from the ATNF Pulsar Database ( Manchester et al.ll2005l ). 
GC MSPs have been excluded because their Ps are usually contaminated by acceleration in 
the cluster potential. Lines of constant characteristic age and surface magnetic field are also 
plotted. 



Table 1. Properties of M22 Used for Scaling 



Ng 5 V esc 


Ly 


logPc 


r c 


(kms -1 ) 


(L Q ,y) 


(L pc- 3 ) 


(PC) 


4 44.7 


2.15 x 10 5 


3.63 


1.2 



Note. - Structural p arameters have been 
taken from lHarrisI (120101 ). and we have used 
D = 3.2 kpc for calculatin g r r in physica l units . 
The escape velocity is from lGnedin et al.l ( 120021 ) . 
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Choose reference GC 

T^-Uv ,ref 






No 


Simulate surveys 







Simulate another 




Yes 



No 



Done 



Fig. 2. — A schematic diagram describing our simulation strategy. The abbreviations for 
model classes are: S = standard models (no scaling); L = luminosity scaling; G = core inter- 
action rate scaling. 
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Table 2. Summary of Models Used to Simulate the Population of NRPs 



Model Class 


N G c 


Scaling 


Sa a 


97 


None 


Sb b 


81 


None 


L c 


156 


Luminosity 


G c 


143 


r c 



Note. - Ten runs were per- 
formed for each value of a v for each 
model class. 

a Sa models use all the GCs stud- 



ied in iPaper I . 



b Sb models exclude GCs with 

L i.4GHz,min > 7 - 5 mJykpc 2 be- 
cause these shallow surveys were 
not very constraining. 

C M22 was used as the reference 
cluster when scaling. 



Table 3. Parameters Used in the Model of Galactic Potential 



Component 


M 


a 


b 


ft 


hi 


P2 


h 2 


03 


h 3 




(10 9 M ) 


(kpc) 


(kpc) 




(kpc) 




(kpc) 




(kpc) 


Disk-Halo 


145 


2.4 


5.5 


0.4 


0.325 


0.5 


0.090 


0.1 


0.125 


Bulge 


9.3 





0.25 


1 

















Nucleus 


10 





1.5 


1 


















Note. — See ^3.21 for parameter definitions. 
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Table 4. Large-area Survey Parameters 



Survey 


£ a 
J sky 


G 


T 

sys 


Vohs 


Av 


^samp 


tint 


P 


Sl374 MHz,min b 




(%) 


(K Jy- 1 ) 


(K) 


(MHz) 


(MHz) 


(fJS) 


(s) 




(MJy) 


PMSURV 


4.6 


0.6 


25 


1374 


288 


250 


2100 


1.2 


160 


P-ALFA 


10.9 


8.5 


25 


1374 


300 


64 


268 


1.1 


100 


GBNCC 


19.2 


2.0 


46 


350 


100 


81.92 


120 


1.1 


2000 


GBTALL 


86.0 


2.0 


46 


350 


100 


81.92 


120 


1.1 


2000 


SKA 


100 


140 


25 


1374 


512 


50 


2100 


1.0 


0.45 



Note. — All surveys use two summed polarizations. 
a Fractional sky coverage 

b Approximte limiting flux density scaled to 1374 MHz assuming a spectral index of —1.6 



L70 




Fig. 3. — The position of all escaped pulsars in Galactic coordinates (as viewed from Earth) 
for model L70. Other models have a similar spatial distribution of pulsars but with varying 
numbers. The Galactic center is at the center of the plot. Globular clusters are represented 
by red diamonds, while pulsars are shown as black dots. 
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Table 5. Simulation Results 



Model TD a 


N b 

-i* esc 


Jv ret 


iv cmit 






Jv det 
















PAT.AFA 




CIRTAT ,T, 

VjJ_l_> _1_ JT. J_J J_J 


SKA 


Sa5 


1.3 x 10 6 


9.5 x 10 B 


1.5 x 10 5 


3 








83 


1.1 x 10 4 


SalO 


1.0 x 10 7 


9.5 x 10 5 


1.2 x 10 6 


40 


6 





2.6 x 10 2 


9.2 x 10 4 


Sa20 


8.0 x 10 7 


9.7 x 10 5 


9.1 x 10 6 


3.9 x 10 2 


51 





3.5 x 10 3 


6.9 x 10 5 


Sa30 


2.7 x 10 s 


1.0 x 10 6 


3.1 x 10 7 


3.0 x 10 2 


7.0 x 10 2 





7.9 x 10 3 


2.4 x 10 6 


Sa50 


1.2 x 10 9 


9.2 x 10 5 


1.4 x 10 s 


5.4 x 10 3 


5.9 x 10 2 





5.4 x 10 4 


1.1 x 10 7 


Sa70 


3.4 x 10 9 


9.7 x 10 5 


4.2 x 10 s 


1.0 x 10 4 


2.4 x 10 3 


1 


1.1 x 10 5 


3.0 x 10 7 


SalOO 


9.9 x 10 9 


9.3 x 10 5 


1.1 x 10 9 


6.1 x 10 4 


5.4 x 10 3 


35 


3.7 x 10 5 


8.3 x 10 7 


Sb5 


3.4 x 10 4 


4.7 x 10 5 


3.5 x 10 3 











2 


3.8 x 10 2 


SblO 


3.8 x 10 5 


4.8 x 10 5 


3.6 x 10 4 


6 


6 





23 


4.3 x 10 3 


Sb20 


3.3 x 10 6 


4.9 x 10 5 


3.3 x 10 5 


63 


62 





2.1 x 10 2 


3.8 x 10 4 


Sb30 


1.1 x 10 7 


5.1 x 10 B 


1.1 x 10 6 


2.5 x 10 2 


1.9 x 10 2 





5.0 x 10 2 


1.3 x 10 5 


Sb50 


5.2 x 10 7 


5.2 x 10 5 


4.9 x 10 6 


9.0 x 10 2 


7.2 x 10 2 





2.4 x 10 3 


5.8 x 10 5 


Sb70 


1.4 x 10 s 


5.6 x 10 5 


1.3 x 10 7 


2.6 x 10 3 


2.8 x 10 3 


2 


7.0 x 10 3 


1.6 x 10 6 


SblOO 


4.1 x 10 8 


5.1 x 10 5 


4.1 x 10 7 


7.6 x 10 3 


6.5 x 10 3 


40 


2.8 x 10 4 


4.4 x 10 6 


L5 


2.0 x 10 2 


1.5 x 10 4 


42 














2 


L10 


8.8 x 10 2 


1.4 x 10 4 


1.7 x 10 2 














10 


L20 


5.0 x 10 3 


1.3 x 10 4 


9.7 x 10 2 





0.17 








58 


L30 


1.6 x 10 4 


1.5 x 10 4 


3.2 x 10 3 


1 








1 


1.9 x 10 2 


L50 


6.9 x 10 4 


2.0 x 10 4 


1.4 x 10 4 


5 


1 





6 


9.4 x 10 2 


L70 


1.9 x 10 5 


2.7 x 10 4 


3.8 x 10 4 


14 


3 





16 


2.7 x 10 3 


L100 


5.9 x 10 5 


3.4 x 10 4 


1.2 x 10 5 


38 


9 


0.57 


47 


7.8 x 10 3 


G5 


53 


7.3 x 10 4 


7 














1 


G10 


1.1 x 10 3 


7.2 x 10 4 


2.0 x 10 2 


0.11 











15 


G20 


1.7 x 10 4 


7.0 x 10 4 


3.3 x 10 3 


3 


0.52 





1 


2.5 x 10 2 


G30 


6.0 x 10 4 


9.7 x 10 4 


1.1 x 10 4 


9 


1 





3 


8.5 x 10 2 


G50 


3.2 x 10 5 


1.9 x 10 B 


6.0 x 10 4 


49 


4 





19 


4.6 x 10 3 


G70 


1.0 x 10 6 


2.6 x 10 5 


1.9 x 10 5 


1.6 x 10 2 


8 





58 


1.4 x 10 4 


G100 


3.3 x 10 6 


3.2 x 10 5 


6.0 x 10 5 


4.6 x 10 2 


24 


1 


2.3 x 10 2 


4.3 x 10 4 



a Model IDs include model class and c v (e.g. Sa5 refers to model class Sa and c v = 5). 

b 7Veac and N ve t are all escaped and retained pulsars produced in the last 200 Myr. These numbers scale linearly 
with the chosen timescale. 

c iV < j et has already been corrected for a 10% beaming fraction. 

Note. — The relative success of some surveys (such as the PMSURV, PALFA, and GBTALL) depends slightly 
on model class. This is because some model classes exclude certain clusters. Pulsars from these clusters will be 
concentrated in a certain region of the sky, but surveys do not have identical sky coverage. 
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Table 6. Birth Rates of GC NRPs 



Model ID 


TZ 


^exhaust 




(psr century" 1 ) 


(Myr) 


Sa5 


1.1 


19 


SalO 


5.5 


3.8 


Sa20 


40 


0.52 


Sa30 


140 


0.15 


Sa50 


600 


0.035 


ba70 


1700 


0.012 


SalOO 


5000 


0.0042 


bb5 


0.25 


83 


SblO 


0.43 


49 


Sb20 


1.9 


11 


Sb30 


5.8 


3.6 


Sb50 


26 


0.8 


Sb70 


70 


0.3 


SblOO 


210 


0.1 


L5 


0.0076 


2800 


L10 


0.0074 


2800 


L20 


0.009 


2300 


L30 


0.016 


1400 


L50 


0.044 


470 


L70 


0.11 


190 


L100 


0.31 


67 


G5 


0.037 


570 


G10 


0.037 


570 


G20 


0.044 


480 


G30 


0.079 


270 


G50 


0.26 


82 


G70 


0.63 


33 


G100 


1.8 


12 



a The time scale to exhaust a supply of 
— 2.4 x 10 5 WDs via ECS, assuming a 
constant birthrate (see text). 



